Predictions from conformal algebra for the 
deeply virtual Compton scattering. 
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Abstract 

Basing on the constraint equalities which arise from the algebra of the collinear conformal 
group and the conformal operator product expansion, we predict the solutions of the leading 
order evolution equations for the non-forward distributions, including transversity, in terms of the 
conformal moments; next-to-leading order flavour singlet coefficient functions for the polarized 
and unpolarized deeply virtual Compton scattering; as well as the contribution from renormalon 
chains to the eigenfunctions of the exclusive non-singlet evolution kernel. 



1 Introduction. 



It is well known that the conformal twist-2 operators possess at leading order (LO) a diagonal 
anomalous dimension matrix 7 = {jjf}, which is, in general, triangular due to Lorentz invariance. 
One may expect that in the conformal limit of the theory (j3 = 0) the former property will 
be fulfilled beyond one-loop order; however, this is not the case and the conformal operators 
mix already in the next-to-leading order (NLO) of perturbation theory. A closer look at the 
conformal Ward identities (CWI) in LO for the flavour non-singlet channel |], 0], which include 
both dilatation and special conformal transformations, shows that the conformal covariance is 
destroyed in the minimal subtraction scheme (MS) by an off-diagonal special conformal anomaly 
matrix 7 C (/) = {7|j'(0}> which depends on the spin I carried by the considered twist-2 operator. 
To restore it to all orders of perturbation theory one is forced to perform a finite renormalization of 
the conformal operators leading to the conformal subtraction (CS) scheme. Assuming a non-trivial 
fixed point g*, i.e. that (3{g*) is zero, this restoration is assured in any order of the perturbation 
theory by conformal constraints which arise from the algebra of the collinear conformal group: 

[D, KJ\_ = %K_ [a(/)+ 7 c (/), 7 ] = 0, (1) 

[P+,K_}_ = 2z(£> + M_ + ) 7 ^ + i)_ 7 ^) = _2 7) ( 2 ) 

where ajf(l) = a(j,l) • Sjf and a(j,l) = 2(j — + 1 + 3). For the case of 0^ theory, this 
constraints can be extended without any assumptions to nonvanishing /5-function 

0, (3) 



a(/)+ 7 c (0 + 2-&(/),7 
9 



where the b matrix is defined as[] 

b jf (l) = 9 jf {2(1 + f + 3)S jf - [1 + (-l)'-'W + 3)} . (4) 

Hence, the contribution to the off-diagonal matrix elements of the anomalous dimension matrix 
due to /5-function is predicted too. Under the conjecture that only gauge invariant composite 
operators contribute to the constraints for the gauge independent anomalous dimension matrix 
of the conformal operators, which turns out to be quite natural, formally the same equation can 
also be derived in Abelian gauge field theory for the flavour non-singlet channel as well as 
for the quark-gluon singlet sector studied in the present paper. The knowledge of the conformal 
anomaly matrix allows to reconstruct the off-diagonal part of the anomalous dimension matrix for 
the corresponding channels. 

1 We use the common definition of the Kronecker symbol Sjji and the following definition of the (discrete) 
f 1, 3-f>0 
k 0, j ~ f < ' 
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Few remarks are in order at this point. Although the above mentioned commutator constraints 
are derived, for the reasons of technical simplification, in the theory with £7(1) symmetry group, 
nevertheless, they are valid in NLO approximation for QQ- and QG-channels we have considered, 
also in the non-Abelian theory. This is true since in the conformal limit, the symmetry breaking 
parts of the kernels do not contain the Casimir operator Ca of the adjoint representation of SU(3). 
Moreover, even in the theory with non-vanishing (3 the correct results can be reconstructed by 
substituting the QED /^-function by the corresponding QCD value [[| and, therefore, this treatment 
is sufficient to derive the reliable NLO predictions in QCD. Thus, some of the results of this 
paper can be regarded as a first step towards the complete derivation of the a s -corrections to 
the exclusive flavour- singlet evolution kernels which are presently completely unknown. Using the 
same technique, the purely gluonic case will be clarified in the nearest future |§. The knowledge 
of all other channels will allow one to check the result with a construction of the (N = 1) SUSY 
Ward identities. 

In this paper we mainly concern the predictions arising from the conformal algebra for the 
deeply virtual Compton scattering (DVCS) process, namely, the solution of the evolution equa- 
tions, which govern the Q 2 - dependence of the corresponding structure functions, and the NLO 
corrections to all singlet coefficient functions in the MS renormalization scheme. The final section 
is devoted to the issue of the renormalon chains contribution to the Efremov-Radyushkin-Brodsky- 
Lepage (ER-BL) kernel which can be incorporated in a straightforward way in our approach. 

2 Evolution of twist-2 conformal moments at LO. 

In the past 0] and more recently 0, §|, different authors introduced non-forward distributions, 
which describe the non-perturbative inputs in certain processes such as DVCS and exclusive 
electroproduction of mesons. In the consequent discussion we accept the following definition of 
the non-forward distribution functions || [7| 

i O r (A, y u) = (ti\<j)*(fin)r<5>{fm,\n]<j) t (\n)\h) = J dxe l ^ x - ~ lXx i O r (x,C), (5) 

in terms of the light-cone Fourier transformation. Here i = Q,G runs over the parton species, 
r corresponds to different Dirac or Lorentz structures, depending on the spin of the constituents 
involved and $ is a path ordered exponential. The momentum fractions of incoming x and outgoing 
x — ( partons are the Fourier conjugated variables of the light-cone positions A and \x. Parameter £ 
is a skewedness of the distribution defined as a +- component of the t-channel momentum A + = £ 
(see next section for discussion of the kinematics of the DVCS process). 

The scale dependence of the light-ray operators, standing on the left hand side of Eq. (|), is 
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governed by the renormalization group equation which is of the generic form (y = 1 — y) 



dy dzK. r (y J z)O r (Xz + fiz,fj,y + Xy). 
Zir Jo Jo 



(6) 



The evolution kernels for the parity-even (V) and parity-odd (A) chiral-even light-ray operators 
are known in the literature |5], ||, [/], ^ ||. Here we supplement this list by the kernel for the 
non-forward transversity (T) distribution [|lO|j : 



K T (y,z) = C F {6(y) 



S(z) 



S(y) - 5(z) + -5(y)5(z) 



(7) 



By exploiting the conversion formula derived in Ref. |Tl| we can easily write the corresponding 
evolution equation in the momentum space 

d_ 

dy 



V 2 J- 2 O r (x,C) = -^Jdx>K r (x,x'X)O r (x'X), 



(8) 



withg 

K T (x,x',C) = C F 



X -^0 / „ /\ i 00 C r*A> 



-Q n (x, x — x') H ]®ii( x ~ C? x — x ') 



x — X 



+ -8{x' - x) 



(9) 



where B^x, x') = (0(x) — 9(x'))/(x — x'). The corresponding ER-BL kernel = 1] 



V T (x,x') = -K T (x,x',l) = C F \\6{x' - x)-—^— + 6(x - x')-—^— 



and the DGLAP kernel (( = and z = x/x') 

P T (z) = -K T (zx', x', 0) = C F | 2 



+ -5(x-x')\,(10) 



follow immediately from the generalized one ®. 

In order to find the solution of the evolution equations of the type given in Eq. (||) we are 
looking for the conformal spin expansion in terms of reduced non-forward expectation values of 
conformal twist-2 operators which are related to the Gegenbauer moments of the corresponding 
distributions, namely, for quark function we have 



C j J dxC} - Q O r (x,C) = (h'fO^h). 



;i2) 



Cj are the Gegenbauer polynomials ]r2| with the normalization fixed below in Eq. ([TTD . The 
gluon case differs only in the change of the index of the orthogonal polynomial by one unity and 



2 The explicit expressions as well as integral and algebraic properties of the O-functions introduced in the main 



text: 91 



^£1,2:2, ...,x n ) = JZ 



da „,m TT n 



Ylk=i ( ax k — 1 + *0) k can be found in Ref. (11 
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shifting j to j — 1. Then, the complete sets of conformal twist- 2 operators for the flavour singlet 
chiral-even un- and polarized as well as chiral-odd channels read: 







M 


< QQ A 




7+75 


QqT 

Ik > 




^+±75, 




D 







^illf (13) 



Here d =d + d and D—d — d —2igT a B a and we suppress the flavour and colour indices for 
simplicity. The operators we have introduced above mix under the renormalization and their 
anomalous dimension matrix 



7 



7 



(14) 



possesses the square blocks lA 7 which are triangular matrices. Since the special conformal anomaly 
does not enter in the conformal constraint ([!]) at LO, the latter takes the simplified form: 



a(7) ' 
o(7) 



'QC 
GC 



QQ. 



1 



GG, 



7, 



0. 



(15) 



The matrix a(l) is diagonal, then it follows from the above equality that all entries in 7 are diagonal 
too and coincide with the moments of the DGLAP evolution kernels^. Therefore, we can easily 
write the solution of the renormalization group equation (§) in terms of the conformal moments, 
for instance, for the quark non-forward distribution in the flavour non-singlet channel: 



NS e^(*,o = £ 



1 



^■(i)c +i c 



- 1 



- cj 



2*-l 



' a.(Q) 
a s {Qo) 



(hT S O r J3 \h) (16) 



with fa = 2JV> 



^tCa and the normalization factor 



-4i/+l 



r 2 (i)r(2z, + J ) 



(17) 



The quark-gluon mixing problem in the singlet channel can be easily solved by straightforward 
diagonalization in the same way as in the case of the deep inelastic scattering (DIS). Note that 
Eq. (|T6|) is valid only for the distributions with a support < x < C 0- 



3 NLO coefficient functions for DVCS. 

In the previous section we have dealt so far with the LO analysis and have seen that the tree 
conformal operators flT5| ) diagonalize the anomalous dimension matrix. As we have mentioned in 

3 For the QG and GQ sectors this is true up to some common factors (see Ref. [^3) and Eq. (|2^) below). 
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the Introduction, this property breaks down in the NLO in the MS scheme even in the theory with 
j3 = 0. However, the conformal constraints (|l|) and (fj) ensure the existence of the CS scheme, 
in which the conformal covariance is restored in the conformal limit of the theory to any order 
of perturbation theory. In this scheme the OPE in terms of the conformally covariant operators 
[ 111 ] holds true in the interacting theory and provides powerful restrictions on the corresponding 
Wilson coefficients and can be served, for instance, as a tool for the calculation of the scattering 
amplitude of two-photon processes in the generalized Bjorken region. To compare the predictions 
with the calculations performed in the MS scheme, the following transformation for the operators 
is necessary 

On = E B Sj ,O jn . (18) 

j'=0 

In the CS scheme the operators Oji in different conformal towers do not mix with each other up 
to terms proportional to (3. The transformation matrix B, determined completely by the special 
conformal anomaly: 

6 = ttV; = 1 - J ^ c + WW - ■ ■ ■ . ( 19 ) 

1 + 3 7 

with J A = 6jjt + iAjjt /a(j, j'), specifies the partial conformal eigenf unctions of the evolution ker- 
nels. This matrix was calculated in Ref. for the non-singlet parity-odd channel and used for 
a consistency check between the available NLO result for the coefficient function of the 7*7 — > 71 
transition form factor and the forward NLO Wilson coefficients for polarized DIS [Dj]. 



Recently, the ct^-correction to the non-forward coefficient function in the singlet channel has 
been calculated for unpolarized DVCS JOJ. In this section we evaluate the one- loop coefficient 
functions for the polarized {T^) and non-polarized (TX\) well: 

1 T F (uj, x, C, Q 2 \a s ) = 'Tfo (u, x) - ^Tfa (u, x, £ Q 2 ) + 0(a 2 s ). 

The predictions for the coefficient function can be obtained by employing the conformal OPE from 
the one side and the conformal constraints from the other. In NLO the symmetry breaking, which 
is proportional to the /3-function, arises only in the off-diagonal part (in the basis of Gegenbauer 
polynomials) of the anomalous dimension matrix 7, not in the coefficient function. So, we are 
able to make the predictions for the NLO coefficient functions without restrictions. 

To overcome difficulties with gauge invariance and Lorentz decomposition of the hadronic 
tensor for the DVCS, we deal, in what follows, only with appropriate projections of the latter: 

^>,C) = iV^J d*ze^(ti\T{J,(0)J u (z)} \h), (20) 
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where = g^ u corresponds to unpolarized scattering and Vt v = ie^ U Q P / (Qp) to polarized case. 
Here l/u = -Q 2 /2(Qp) and Q = q - A. 

On the one hand the amplitude for DVCS can be represented in the factorized form as convo- 
lution of the non-forward distribution ([5]) and the perturbative coefficient function 

F r (uX)= fdx J2 l T r (u,xXiQ 2 \a s yO r (xX). (21) 

i=Q,G 

For the hand-bag diagram^ we have to LO in the coupling ^T^(u, x) = uj/{xuj — 1). 

On the other hand the conformal OPE |TJ], |15| for the product of two electromagnetic currents 
sandwiched between the hadronic states, in the kinematics appropriate to the process in question, 
gives the following prediction: 



r r { u , o = E E 'qf K Q 2 /^ 2 ) " j+l 2Fi 

i j=Q 



1+j + i-yf, 

4 + 2j + 2 <7f 



uC)(h'\%\h). (22) 



Here 2 -Pi is confluent hypergeometric function [17|, the index i sums up the eigenvectors l O F jj of 



the singlet channel, which possess the diagonal anomalous dimensions , and 8 Cj is the Wilson 
coefficient of the forward scattering. We have the convention that in the interacting theory the 
scale dimension is shifted by the amount 2 l/ y F compared to canonical one: dj = d^ n + 2 l/ y F . In Eq. 
(p2|) fi 2 is a renormalization scale, while we have put the factorization scale [i 2 F equal to Q 2 . Notice 
that for practical purposes MS schemes will be used where the manifestly conformaly covariant 
form of the OPE ([22]) is hidden. Fortunately, the product of two electromagnetic currents is a 
renormalization group invariant, so that Eq. ( plf ) computed in any scheme is equivalent to Eq. 
22|) derived in the CS scheme. 



Since the a s -corrections in Eq. (21) come from the NLO result for the coefficient function and 
eigenf unctions of the evolution kernels we are able to predict the former from the combined use of 
the conformal OPE ( p2"D and the special conformal anomaly which fixes precisely the corrections 
to the latter. The general conformal decomposition of the non-forward distributions looks like 

00 

i O(x,()= E T,%M\as){h'\ k O jj \h), (23) 

k=Q,Gj=0 

where the partial conformal waves are generalized, beyond one-loop level, to non-polynomial 
functions which are the subject of the constraints. For the quark case, which is required for the 
present application, they have the form 

Qk <PMX\a s ) = [5 Qk 5(x -y) + %* k *(x,v,Q} ® J^n\ f 1 " fj q (4 - l) . (24) 



4 We omit everywhere the crossed u-channel contribution. 
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where v = i/(k) and we introduced the shorthand notation ® = / dy. Here w $ is known from 
Ref. 0, while Q G (p can be found with the help of the special conformal anomaly matrix in the 
quark-gluon channel. Fortunately, in this sector no subtle conformal symmetry breaking effects 
occur as it happens in the quark-quark transition amplitude. Thus, the special conformal anomaly 
can be obtained immediately and reads: 

QG^c) = _^S(Z) QG 7 . (25) 
2ir 

Now we have all necessary results to solve the problem in question. First we outline the proce- 
dure for the evaluation of the coefficient function for the polarized non-forward gluon contribution 
and then we just list the results for other cases since the machinery we have developed is quite 
general and can be easily extended to the quark coefficient functions as well. 

3.1 Polarized case. 

For the QG-sector we find from Eq. (|25|) that the correction to the eigenf unctions are completely 
expressed in terms of the shift operator: 

QG <P(x, y, C) = (X - V)S(x, z) ® Q G K A (z, y, (), (26) 

which generates the shift of the Gegenbauer polynomials index: 

[xx]^ +p C" +p (2x-l). (27) 



S{x,y)®[nr-tC]{2y-l)=± 



p=0 



In Eq. (|26|) X is an identity operator and D extracts the diagonal part in the expansion of any 
test function r(x,y) with respect to Cj, i.e. 

oo 

Vr(x,y) = Y\xx] v ^/N j {v)q{2x - l)r jj q(2y - 1). 

The generalized evolution kernel^ 

QG K A (x, x', C) = 2N f T F Q\ 12 (x, x-(,x-x') (28) 
can be diagonalized with the following identity 

C / dyCj (2 V - - l) <*^(y, x, C) = QG ltcU ^ - l) with = - ^J F + (29) 



5 These generalized kernels posses the following limits: for ER-BL kinematics x < £ : K(x,x',() 
-V{x/(,x' /()/(, for froward DGLAP evolution C = 0: K{x,x',0) = -P{x/x'). 



7 



The coefficient ^ G ^ A in front of the Gegenbauer polynomial coincides with the DGLAP moments 
of the kernel ® G K (y, x, 0) up to the common factor 6/j, which arises as a result of the conventional 

3 

definition of the Gegenbauer polynomials, namely, the coefficient of x J in CHx) is 3/j times the 

5 

coefficient of in C^_ 1 (x); an additional factor of 2 comes from the argument of the polynomial. 

For the off-diagonal part of the amplitude we get from the conformal expansion fl2"2] ) the 
following expression 

(J - V) [ Q T (0) (w, y) ® S(y, z) - Q T {0) (u, z) ln(l - zu)] ® QG K A (z, x, () ® g O a (x, () (30) 



3=0 



d 
dp 



l F 1 



p=0 



+ p, 1+ 3 + p 
A + 23 + 2p 



ujn{h'\ G o A 3] \h). 



Its normalization is related to the LO prediction of the conformal OPE for the DVCS scattering: 

'1+3, 2+3 



F A Q {u, C) = -2 E B (J + U + 2)^' +1 2 Fi 



3=0 



A + 2j 



uJ 



C)(h'\ ( ^0 A jJ \h), 



(31) 



To fix the diagonal part of the hard scattering amplitude we have to put = in Eq. (|22"D, 
but keep the NLO result for the forward coefficient function of the polarized gluon distribution 
§ G C/(g 2 //i 2 ) 0, |TJ|. We can write the result convolution 



Q T m (u,z)®F(z,x,()® G A (x,() 

oo 

= -2j2 G C A (Q 2 /p 2 )B(j + 1,3 + 2V +1 2 Fi 

3=0 



(32) 



1+3, 2 + j 
4 + 2j 



uj()(h>\ G A n \h), 



with the following explicit expression for the function F, which corresponds to the DIS coefficient 
function in the MS scheme [ [19|j : 



6 



Q 2 x 



2N f T f - \ (2x - 1) In + 3 - Ax \ <g> x 



prx 



(33) 



Q 



F(z,x,C) = - QG K A (z,x,C)^-+r + [1 -T?ln(l - zu)\ QG K A (z, x, (). 



/' 

By adding the results of Eqs. (|30|) and (|33|) together and subtracting the solution (|26|) of the 
evolution equation for (3 = 0, we finally come to the coefficient function in the MS scheme: 

V2 



G T A (u,x,, 



C, Q 2 ) = Q T {0) (u, y) ® [ QG K A (y, x, () In *L + [l n (1 - yu ) - 1] Q G K A (y, x, C)}. (34) 



As an outcome of the same but a little bit more involved analysis for the quark scattering 
amplitude we get: 



®T A (u, x, (,Q 2 ) = %) (u,y)®\ ^K A (y, x, () In % - \ \ Q K b (y, x, ()] , + |<% - x) 



m ~ I ' ^ - p? 2 

+ hi(l-yu;)^K A (y,x,0-\ Q G(y,x,0 



(35) 
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Here 



QQ K A (x,x',() = K T {x,x',() + C F 0° ul {x,x-(,x-x'), 



Q G(x,x' } () = ^—^\n 



x ,_ ( - .e u {x-{,x-x') + ^ _ ; 



(36) 

e^^x -x')(37) 



(38) 



and a new kernel ®G(x, x', () comes from the expression for Q^cp, which is of the form for (3 = 0: 



x, y, C) = (l-V)\ S(x, z) ® ^K A (z, y, () - \^G(y, x, () 



(39) 



In the limit to = 1 and ( = 1 we obtain from Eq. ( pq ) the results of Refs. [20, 21, 22] for the pion 
transition form factor F 1 



7T77* 



3.2 Unpolarized case. 

Since the QQ-evolution kernels are the same in the polarized and non-polarized cases, the only 
difference will arise from the forward coefficient functions [p3|| . Straightforward calculation gives 
for the quark as well as for gluon non-polarized coefficient functions: 

Q 2 5 u 



Q T& («, x, C, Q 2 ) = Q T (0) (u,y)®\ QQ K v (y, x, () ln ^ - - [^K b (y, x, ()} + + 5(y- x) 



+ [ln (1 - yu) - 1] QQ K v (y, x, C) - [ Q G(y, x, () 

Q 2 1 r, 

+ ln(l-yu)Q G K v (y,x,()\, 



(40) 



G T ( ^, a;, C, Q 2 ) = Q T {0) (uj, y) ® \ QG K v (y, x, C) In ^ - ~ [ QG K v (y, x, C) + QG K A (y, x, C) 



the momentum space evolution kernels °-(Q> G )K v read 
^(x,x',C) = QQ ^ A (x,x',C), 

QG ir y (x,x',C) = QG K A (x, x', C) - ANfTp ^'^ e° ul (x, x-Cx-x 1 ). 



(41) 



(42) 



These expressions coincide with the ones obtained in Ref. |T6| provided we have made some 
obvious redefinition of their conventions. Taking the limit intrinsic to the exclusive kinematics we 
obtain the well know result [20|, for the radiative corrections to the coefficient functions of the 
parity-even vector meson transition form factor. 
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4 Renormalon chains in the non-singlet evolution kernel. 



The last issue of this paper that can be studied by the presently developed methods is the renor- 
malon chains resummation of the ER-BL evolution kernel. Recently, this problem has been con- 
sidered in the MS scheme for the non-singlet channel (see also the work for the DGLAP 
splitting kernels |26|, |27|). It has been found that the resulting series possesses a nonzero conver- 
gent radius, so that the infrared renormalons are absent in the kernels. From our point of view this 
result can be achieved by a partial resummation of conformal anomalies. The 5-matrix defined 
in Eq. ( |T9| ) diagonalizes the anomalous dimension matrix for the case (3 = 0. From the confor- 
mal constraints ([3]) we know that in the MS scheme its off-diagonal part is induced by 7 C + 2^6. 
Therefore, the operator 

B = t~ (43) 



9 

diagonalize the anomalous dimension matrix 6 . This formula can be understood as a resummation 
prescription for the conformal anomalies that induce the off-diagonal part. Knowing this matrix 
one can construct the eigenfunctions (f>i{x\a s ) = QQ 0;(x, l\a s ) of the ER-BL kernel 

00 xx £ 

UAocs) = E -W7sS C I( 2x - !)%(«.)■ ( 44 ) 
j=i iy i\2> 

In Ref. a differential equation for the eigenfunctions was derived in the conformal limit. Re- 
peating the steps, it is easy to restore the neglected (3 term. Employing the relation (see Appendix 
in |J) 

a(j, l)B 3l (a s ) = - | (^(0 + 2^6(Z)) B j (45) 

and the eigenvalue equation for the Gegenbauer polynomials [|l^| we are able to derive an integro- 
differential equation for the eigenfunctions in the full theory, valid for the MS scheme: 

xx-^4>i(x\a s ) + (I + 1)(/ + 2)<pi(x\a s ) = - dy ijf(x,y,a s ) + 2-bi(x,y)j^i(y\a s ), (46) 

where the kernels can be obtained from their non-forward analogues introduced in the preceding 
sections and read 

^(x,y,a s ) = ^{2j 1 dzb l (x,z)[WV(z,y)) + + [g(x,y)] + } + O(a 2 s ), (47) 



6 Since the B-matrix depends on the running coupling, the diagonalization of the anomalous dimension matrix 
provide the solution of the renormalization group equation only in the fixed coupling regime. 
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k(x,y) 
} V(x,y) 



2[l + l 



2x - 1 d 



C F 9{y - x) 



2 dx 
ln(l - 



8{x - y), 



v 1 



y-x 



C F 6(y 



*)- 1 + - 

y V v 



+ 
i 



J a; — >• 


X 




y 




X 




[y 



X 

y . 



(48) 
(49) 
(50) 



The non-perturbative solution of Eq. fl46|) is unknown. However, in the large Nf limit only the 
term proportional to (3 survives. Taking into account also the non-leading Nf terms included in 



the /3-function, the equation (|4q) can be reduced to 

A 2 

dx 2 



xx— (j>i{x\a s ) + {l + I) (I + 2) 



Mx \a s 



(51) 



- [ dy 
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211 + 1 + 



2x- 1 d 



dx 



5(x -y) + 0{a s ) \ <f>i{y\a s 



Taking the known LO contribution on the right hand side yields a differential equation which has 
the solution: 



(pj(x\a s ) ~ [xx]*~b C- 9 (2x - 1). 



(52) 



This result coincide in the large Nf limit with the resummation of the vacuum polarization bub- 
bles in Refs. [[24], |25| and verifies the 'naive non-abelianization' hypothesis employed there |25| . 
Moreover, it is very interesting to see that this hypothesis can be extended in our case to the whole 
/3-function. Unfortunately, the eigenfunctions (|52D have no practical importance for the solution 
of the evolution equation. Obviously, in this treatment the special anomaly matrix contribution 



is completely discarded, however, it was demonstrated by numerical studies in Ref. |f28| that a 
considerable cancellation occurs between the former and the symmetry breaking term due to the 
/3-f unction for the real case Nf = 3. 



5 Summary. 

In the present study we employ the powerful restrictions from the conformal algebra on the 
amplitudes in the massless theory. In the leading order of perturbation theory it enables us to 
give the solution of the evolution equations for the off-forward distribution functions in terms 
of the conformal partial wave expansion. Beyond the leading order the conformal symmetry 
breaking term induced by nonvanishing /5-function does not appear in the coefficient function so 
that we can rely on the hypothetical conformal limit of the theory in order to make the predictions 
for the two-photon processes in the generalized Bjorken kinematics. To this end the knowledge 
of the forward Wilson coefficients in the OPE and the diagonal anomalous dimensions of the 
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corresponding operators which were determined earlier from the study of deep inelastic scattering 
is required . 

The development of the present method is far from being finished since up to now all the results 
were derived in the theory with 2/(1) symmetry group. The extension of the conformal constraints 
to purely Yang-Mills theory M will allow to get much deeper insight into the conformal symmetry 
breaking by the non-abelian fields. 
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ported by the Russian Foundation for Fundamental Research, grant N 96-02-17631 and Deutsche 
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